Exact dissipative cosmologies with stiff fluid 
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^ ! I. INTRODUCTION 

o: n 

Ps| i Over thirty years ago Misner M suggested that the observed large scale isotropy of the Universe is due to the 

,— H ' action of the neutrino viscosity which was effective when the Universe was about one second old. Dissipative bulk 
^ , viscous type thermodynamical processes are supposed to play a crucial role in the dynamics and evolution of the early 
^^ ' Universe. There are many processes capable of producing bulk viscous stresses in the early cosmological fluid like 
Cn ] interaction between matter and radiation, quark and gluon plasma viscosity, different components of dark matter, 

^ • particle creation, strings and topological defects etc. [Q. 

Traditionally, for the description of these phenomena the theories of Eckart |S] and Landau and Lifshitz B were 

used. But Israel |^ and Israel and Stewart |6| have shown that the Eckart-type theories suffer from serious drawbacks 

^) I concerning causality and stability. Regardless of the choice of equation of state, all equilibrium states in these theories 

f^ , are unstable and in addition signals may be propagated through the fluid at velocities exceeding the speed of light. 

C ■ These problems arise due to the first order nature of the theory i.e. it considers only first-order deviations from the 

^^ , equilibrium. The neglected second-order terms are necessary to prevent non-causal and unstable behavior. 

J^ ■ A relativistic consistent second-order theory was found by Israel [|| and developed into what is called transient 

^^ \ or extended irreversible thermodynamics 0], H, |p|. Therefore, the best currently available theory for analyzing 

O . dissipative processes in the Universe is the full Israel-Stewart-Hiscock causal thermodynamics. Exact general solutions 

^ '■ of the field equations for flat homogeneous Universes filled with a full causal viscous fiuid source for a power-law 

?H - dependence of the bulk viscosity coefficient on the the energy density have been obtained recently in 0- |ll|. 

^*^' It is the purpose of this Letter to obtain the general solution of the gravitational field equations for the case of 

a bulk viscous cosmological fiuid, with the bulk viscosity coefficient proportional to the energy density, obeying the 

Zeldovich (stiff) equation of state. In this case the solution can be presented in an exact parametric form. The 

behavior of the scale factor, deceleration parameter, viscous pressure, viscous pressure-thermodynamic pressure ratio, 

jrt ■ comoving entropy and Ricci and Kretschmann invariants is also considered. 

II. GEOMETRY, FIELD EQUATIONS AND CONSEQUENCES 

For a Friedmann-Robertson- Walker (FRW) Universe with a line element ds^ = dt^ — a^{t) [dx^ + dy^ + dz^) filled 
with a bulk viscous cosmological fluid the energy-momentum tensor is given by 

Tt^ip + p + n)u,u''-(j> + U)St (1) 

where p is the energy density, p the thermodynamic pressure, 11 the bulk viscous pressure and Ui the four velocity 
satisfying the condition UiU^ = 1. We use units so that SirG = c — I. 

The gravitational field equations together with the continuity equation T*^. — imply 
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3i72 =p,2H + 3H^ = -p-Il,p + 3{p + p)H = -SiJH, (2) 

where H = d/a is the Hubble parameter. The causal evolution equation for the bulk viscous pressure is given by |12| ] 

rU + U=-3^H-^TuUH+^-^-^y (3) 

where T is the temperature, ^ the bulk viscosity coefficient and t the relaxation time. Eq. M) arises as the simplest 
way (linear in H) to satisfy the H theorem ( i.e., for the entropy production to be non- negative, S\ = iP /£^T > 0, 

where 5"* = eN'^ — WT^'^ ^^ ^^^ entropy flow vector, e is the entropy per particle and N^ — nu^ is the particle flow 
vector) §, §. 

In order to close the system of equations @) we have to give the equation of state for p and specify T, ^ and r. As 
usual, we assume the following phenomenological laws |l3]: 

P={l-l)p,i = ap',T = Pp\t = ^p-' = ap'-\ (4) 

where 1 < 7 < 2 , and a>0, /3>0, r>0 and s > are constants. Eqs. (||) are standard in cosmological models 
whereas the equation for t is a simple procedure to ensure that the speed of viscous pulses does not exceed the speed 
of hght. 

The requirement that the entropy is a state function imposes in the present model the constraint r = (7 — 1) /7 PI 
so that < r < 1/2 for 1 < 7 < 2. 

The growth of the total comoving entropy E over a proper time interval {to,t) is given by |13]: 



S](i) - E {to) = -3kg^ / UHa^T^^ 

Jto 



dt, (5) 



where ks is the Boltzmanns constant. 

The Israel-Stcwart-Hiscock theory is derived under the assumption that the thermodynamical state of the fluid is 
close to equilibrium, that is the non-equilibrium bulk viscous pressure should be small when compared to the local 
equilibrium pressure |n| « p = {-j — 1) p |1^]. If this condition is violated then one is effectively assuming that 
the linear theory holds also in the nonlinear regime far from equilibrium. For a fluid description of the matter , the 
condition ought to be satisfied. 

To see if a cosmological model inflates or not it is convenient to introduce the deceleration parameter q = ^^ 1 = 

p+3p+3 ^ rpj^g positive sign of the deceleration parameter corresponds to standard decelerating models whereas the 
negative sign indicates inflation. 

With these assumptions the evolution equation for flat homogeneous causal bulk viscous cosmological models is 

H + 3HH + S^-'a-^H^-^'H - (1 + r)H-^H^ + -h-2)H^ + ^-^jH"^-^' = 0. (6) 

4 2a 

By introducing a set of non-dimensional variables h and by means of the transformations H — a^h^t = ^^9, 

with ao = i^T') ^~^° i * 7^ 1 ^'^"^ using the expression of r as a function of 7, Eq. (g) takes the form 

By denoting n = [1 — 2s) / (1 — r) and changing the variables according to ft, = j/^/(^^'"\ 77 = J y^/ "^ dO Eq. ffl) 
becomes: 

+ (2 + y«)| + (2.-l)j/ + 2;"+i = O. (8) 



III. GENERAL SOLUTION FOR A STIFF COSMOLOGICAL FLUID 



We consider the case of a stiff cosmological fluid with equihbrium pressure equal to the energy density by taking 
the values of the parameters 7 = 2 and r — 1/2. By introducing the substitutions v — 1/u and u — dy/drj , Eq. (@) 
can be transformed to a second type Abel first order differential equation: 



dv 
dy 



(2 + 2/") v^ + y-^+^V 



(9) 



A particular solution of Eq. (O) for rt = 1, s = 1/4 has been obtained in |l4| . We consider now another case of a stiff 
cosmological fluid with bulk viscosity coefficient ^ linearly proportional to the energy density p,S, = ap. Consequently 
s = 1 and n = — 2 and hence Eq. (pi) becomes 



dv 



1 



dy y 



B{y) 



d A{y) 
dy B{y) 



(10) 



where A{y) = 2j/^ — 1 and B{y) = —y. By introducing a new variable a — - ~ ■wt), Eq. (H) can be written in the 
general form 



^ = aB{y)+A{y), 



(11) 



dy ^ 2 1 

-— = 2y -cry-1. 
da 



(12) 



Hence we have transformed the initial Abel type equation into a Riccati differential equation. A particular solution 
of Eq. (12)is given by 



y 



l + a2 



and therefore the general solution of Eq. (IT^ is 

y{a) = -crA(o-) 



A2(a)e-V 



C-2jA^{a)e'^da 



(13) 



(14) 



where A (a) = (l + cr^) and C is a constant of integration. 

Hence the general solution of the gravitational field equations for a Zeldovich causal bulk viscous fluid filled flat 
FRW Universe, with bulk viscosity coefficient proportional to the energy density, can be obtained in the following 
exact parametric form, with a taken as parameter: 



t{a)^to = -—- 



da 



3ao J 2/(ct) 



H (a) = aoy^ (a) , a{a) = ao exp 



-3 / y{<^)dcr 



,p{a) ^ialy^ia), 



q{a) = 5 



3(7 



y2(cr) y(a) 

E (a) = S (ao) + 



,Ilia)^6aly'ia)[y^ia)-ayia)-l], 



j^"T° I y (o-) [y^ (c^) - f^y (c^) - i] exp 



y2 (cr) y (a) 
2 I y (<t) d<j da. 



(15) 
(16) 
(17) 



where ao, to and E (ctq) are constant of integration and ao = 2/3a. 

The singular or non-singular character of the solution for all times i > can be checked from the finite (infinite) 
character of the Ricci invariant BijR^^ and Kretschmann scalar RijkiR''-''^'' , given by 



/ = R^^R^i ^ Ua^y^ [(3 + Say - 5y^) (3 + Say - Ay^) + y*] 



J = RrjkiR''''' = 12a^2/4 (3 + Zay - hy 



„2^2^y4 



(18) 
(19) 



IV. DISCUSSIONS AND FINAL REMARKS 

The evolution of the causal bulk viscous Zeldovich fluid filled flat Universe starts generally its evolution from a 
singular state, as can be seen from the singular behavior of the invariants / and J, presented, for different values of 
the integration constant C, in Fig.l. Generally, the dynamics of the Universe depends on the numerical values of C 
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FIG. 1. Time evolution of the invariants / = RtjR'-' (solid curves) and J — RijkiR^''^^ (dashed curves) for different values of 
the integration constant C (ao = 3/2). 

The behavior of the scale factor is presented in Fig. 2, for some speciflc values of the integration constant. The 
evolution is expansionary for all times. At the initial moment the scale factor is zero, while the energy density tends 
to infinity. 
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FIG. 2. Time variation of the scale factor a for different values of the integration constant C: C = —1.2 (solid curve) 
C — —1.1 (dotted curve) and C — —1 (dashed curve) (qq = 3/2). 



The dynamics of the deceleration parameter, shown in Fig. 3, indicates, for the chosen range of the integration 
constant, a non-inflationary behavior for all times, with g > 0,Vt > 0. 
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FIG. 3. Evolution of the deceleration parameter q for different values of the integration constant C: C = —1.2 (solid curve), 
C — — l.I (dotted curve) and C — —1 (dashed curve) (qq = 3/2). 

The bulk viscous pressure 11 is negative during the cosmological evolution, 11 < 0, Vf > 0, as expected from a 
thermodynamic point of view. In the large time limit, as can be seen from Fig. 4, the viscous pressure tends to zero. 
In the same limit the bulk viscosity coefficient also becomes negligible small. 
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FIG. 4. Time behavior of the bulk viscous pressure IT for different values of the integration constant C: C = —1.2 (solid 
curve), C — —l.I (dotted curve) and C = —1 (dashed curve) (qq — 3/2). 



The ratio I of the bulk and thermodynamic pressures, I = \Tl/p\ is presented in Fig. 5. For all times the condition 
I < 1 holds and therefore the model is thermodynamically consistent. 




t 

FIG. 5. Variation with time of the ratio I = — for different values of the integration constant C: C = —1.2 (sohd curve), 
C — —1.1 (dotted curve) and C — —1 (dashed curve) (qq = 3/2). 

During the cosniological evolution a large amount of conioving entropy is produced, with the entropy E increasing 
in time and tending in the large time limit to a constant value. 
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FIG. 6. Dynamics of the comoving entropy E for different values of the integration constant C: C 
C = —1.1 (dotted curve) and C 



1 (dashed curve). We have used the normalizations qq — 3/2 and — ^ 'a " 



1.2 (solid curve) 
1. 



Due to the specific equation of state satisfied by the bulk viscosity coefiicient, the relaxation time t is a constant 
in the present model. 

Most of the known exact solutions of the gravitational field equations with a causal bulk viscous cosmological 
fluid do not satisfy the required conditions of the thermodynamic consistency, leading to an inflationary behavior 
and violating the condition of smallness of bulk viscous pressure. The solution represented by Eqs. (p^)-(p^) can 
consistently describe the early dynamics of a superdense post-inflationary era in the evolution of the matter dominated 
Universe, when, as expected, the bulk viscous dissipative effects play an important role. 
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